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Abstract 

We give a derivation of the thermodynamic restrictions on the constitutive 
relations of an electrically polarizable and finitely deformable heat conducting 
elastic continuum, interacting with the electric field. This is made following 
the method of Coleman-NoU in a thermodynamic theory with the Clausius- 
Duhem inequality. 

1 Introduction 

We follow the theoretic path and some notations written in [T], where the 
thermodynamics restrictions for an elastic body are deduced in a theory 
where the method of Coleman-Noll with the Clausius-Duhem inequality are 
used and where, differently from here, the heat flux is an independent variable 
that obeys a generahzed Cattaneo's equation. 

This is made in order to set up the first step for extending [1] from 
thermoelasticity to thermo-electroelasticity within a second-sound theory. 
As a particular case we obtain formulae of Tiersten's fundamental paper [2], 
where, to describe the heat conducting continuum in interaction with the 
electric field, the theory is based on a macroscopic model. 



2 Preliminary Definitions 

Let £ denote a three-dimensional Euclidean point space. We consider a 
body B whose particles are identified with the positions X G £^ they occupy 
in a given reference configuration B. 

A positive mass measure is assigned to i? by a referential mass density 
Pr{.) : B ^ (0, oo), so that m(P) := Jp pndV is the mass of the part P of 
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B. 

The material filling B is characterized by a given process class iP(-B) of 
S as a set of ordered 10— tuples of functions on B x 1? 

p = (x(.), e{.), ^(.), e{.), t(.), P(.), q(-), b(.), r(.)) G 1P{B) (1) 

defined with respect to B , satisfying the balance laws of linear momentum, 
moment of momentum, energy, the entropy inequality and the field equations 

of electrostatics, where 

• X = x(X, t) is the motion, 

• 9 — 9(K, t) e (0, oo) is the absolute temperature, 

• ip = (/?(X, t) is the electric potential, 

• e = £(X, t) is the specific internal energy per unit mass, 

• T) — r](K, t) is the specific entropy per unit mass, 

• T — t(X, t) (S = S(X, t)^ is the Cauchy {first Piola-Kirchhoff 
stress tensor), 

• P = P(X, t) {w — IP(X, t)^ is the spatial (referential) polariza- 
tion vector, 

• q = (Q — Q(X) ^)) is the spatial (referential) heat flux 
vector, 

• h — b(X, t) is the external specific body force per unit mass, 

• r = r(X, t) is the radiating heating per unit mass. 

Any motion x(., .) of -B is a regular function, e.g. two-times continuously 
differentiable with respect to t and at each time t continuous and invertible 
from B into £. 

We use Grad, grad, Div, div, to denote gradient and divergence with 
respect to X and x, respectively, whereas a superposed dot denotes the ma- 
terial time derivative. 

The deformation gradient F at X at time t is given by 

F = F(X, t) = Grad x(X, t) , (2) 
and the invertibility of the deformation is sured by the condition 

J = detF > . 
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The velocity v of X at time t is given by 

V = v(X, t) = i(X, t) . (3) 

The law of conservation of mass is expressed by 

pji = pj , p + pdiv\ = , (4) 

where p = p(X, t) is the mass density of X at time t. 

The spatial heat flux vector q and the spatial polarization vector P are 
related with their referential counterparts by 

IP = JF-^P, Q = JF-^q. (5) 

The electric potential if, together with the polarization vector field P, 
determines the eulerian electric displacement field, in Gaussian units, by the 
equality 

D = E*^ + 47rP , (6) 

where E*^ = — V^v^ is the (Maxwellian) spatial electric vector field. 

Note that any two corresponding referential and spatial 'energy-fiux' vec- 
tors, related by 

H = JF-^h, (7) 
have spatial and referential divergences that are related by 

DIVm = Jdivh . (8) 

Hence the referential electric displacement field is 

A = JF-^D = JF-^E*^ + 47rIP . (9) 

The spatial and referential polarization vectors per unit volume are re- 
spectively defined by 

7r = P/p, U = lP/pn. (10) 



3 Spatial Description 

3.1 Local balance laws in spatial form 

Under suitable assumptions of regularity, and using (jl]), the usual integral 
forms of the balance laws of linear momentum, moment of momentum, en- 
ergy, the field equations of electrostatics, and the entropy inequality are 
equivalent to the spatial field equations 

pv = dit;T + P • V^E^^ + pb , (11) 
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pe = T ■ Vv — (izf q + E • p-k + pr , 



(12) 



E 



M 



(13) 



pi] > p{r/6) — div{q/9) . 
Incidentally, note that by the continuity equation we find 



(14) 



■pir = E*' • (P + Pdw^r) . 
Let ijj = ipO be the specific free energy per unit mass defined by 



(15) 



(16) 



Then (1121) and (HM yield the dissipation inequality 



p(V' + r/^) - T • Vv + -q • g + pTT ■ E^ < 



(17) 



where g = grad6(K, t) is the spatial temperature gradient. 

We note that Eqs. ( fTTl) . (|T2l) . ( fT^ and (fT6l) respectively coincide with 
Eqs. (3.23), (3.40), (3.43) and (4.2) of [2]. 

3.2 Constitutive Assumptions in Spatial Form 

Let V be an open, simply connected domain consisting of 4— tuples (F, ^, E^''^, g), 
and assume that if (F, 9, E^^ g) G V, then (F, 9, E^^ 0) G V. 

Assumption 3.1 For every p G JP{B) the specific free energy ip(X., t), the 
specific entropy ri(X., t), the Cauchy stress tensor T(X, t), the specific polar- 
ization vector P(X, t), and the heat flux q(X, t) are given by continuously 
differentiahle functions on T> such that 



V' = V(F, ^,E^^ g). 



(18) 



?7 = 7J(F, e,E^^ g), 
t = t(F, e,E^ g), 
P = P(F, ^,E*^ g), 
q=q(F, ^,E*^g). 



(19) 
(20) 
(21) 
(22) 
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Of course, once p(.), P(.)) ^(0 and //(.) are known, then equality ffTUj) 
gives the continuously differentiable function e{.) determining £:(X, t) such 
that 

e = e{F, e, E^^ g) . (23) 

Also note that the dependence upon X is not written only for brevity; 
when the body is not materially homogeneous it becomes active. 

3.3 Coleman-Noll Method and Thermodynamic Re- 
strictions 

Given any motion x(X, t), temperature field ^(X, t) and electric poten- 
tial field v'(X, t), the constitutive equations (fT8l) -( l2Tl) determine £(X, t), 
r/(X, t), t(X, t), P(X, t), q(X, t), and the local laws and ([12]) deter- 
mine b(X, t) and r(X, t). Hence for any given motion, temperature field 
and electric potential field a unique process p is constructed. 

The method of Coleman-Noll [3] is based on the postulate that every 
process p so constructed belongs to the process class JP{B) of B, that is, on 
the assumption that the constitutive assumptions ( |T8l) -(l22l) are compatible 
with thermodynamics, in the sense of the following 

Dissipation Principle For any given motion, temperature field and 
electric potential field, the process p constructed from the constitutive equa- 
tions ^B^-^2B^ belongs to the process class IP{B) of B. Therefore the con- 
stitutive functions ( I7g|) -( f^) are compatible with the second law of thermody- 
namics in the sense that they satisfy the dissipation inequality ( fyTj ). 

It is a matter of routine to extend to thermo-electroelasticity the remark 
for thermoelasticity written in [3], on page 1119, lines 8-30 from top; here we 
write such extension by paraphrasing Coleman's remark. 

Remark 3.1 Let A(t) be any time- dependent invertible tensor, a{t) any 
time- dependent positive scalar, a(t) any time- dependent vector, P(t) be any 
time- dependent scalar, h{t) any time- dependent vector, and Y any material 
point of B whose spatial position in the reference configuration B is Y. We 
can always construct at least one admissible electro -thermodynamic process 
in B such that 

F(X, t), ^(X, t), g(X, t), E^^(X, t) 
have, respectively, the values A(t), a(t), a, b X = Y. 
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An example of such a process is the one determined by the following de- 
formatoin function, temperature distribution and electric potential: 

X = x(X, t)^Y + A{t) [X - Y] , (24) 

= ^(X, t) = a{t) + [A'^{t)a{t)] • [X - Y] , (25) 

if = ^(X, t) = Pit) + [A^(t)b(t)] • [X - Y] . (26) 

Thus, at a given time t, we can arbitrarily specify not only F, 9, g and 
but also their time derivatives F, 9, g and at a point Y and be sure 
that there exists at least one electro-thermodynamic process corresponding to 
this choice. 

The next theorem can be proved by using this remark. 

Theorem 3.1 The Dissipation Principle is satisfied if and only if the fol- 
lowing conditions hold: 

(i) the free energy response function ip{F, 9, E-^, g) is independent of the 
temperature gradient g and determines the entropy, the first Piola-Kirchhoff 
stress, and the polarization vector through the relations 

r]{F,9,E'') = -de^{F,9,E''), (27) 



t(F, 9, E^) = pF^V'(F, 9, E^) , (28) 

7f (F, 9, E^) = -d^M^iF, 9, E^) . (29) 
(a) the reduced dissipation inequality (Fourier inequality) 

q ■ g < (30) 

is satisfied. 

Proof . By the chain rule we have 

7P^dF^-F + de^-9 + B^m'^ -F + dg^-g. (31) 

Note that 



dv' dv' dX^ dX^ d dx 



(32) 



dxj dX^ dx^ dx^ dt dX^ ' 
that is, 

Vv = F-'^F , (33) 
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and thus 

r.Vv = F-V-F. (34) 

Thus by substituting Eqn. (13T|) together with the constitutive equations 
(|T8l)-( l22|) into the dissipation inequahty (fTTI) . gives 



{pd^iJ - F^^r) • F + {pdetp + r])9 + {pd^Milj + P) • (35) 

+pd^il) • g + iq • g < . (36) 

Now we proceed in the Coleman- Mizel [3| method: by Remark 13.11 we 
can we can state that F, 6, E*^ and g can be assigned arbitrary values 
independently from the other variables and the theorem is proved. <) 

A consequence of the reduced dissipation inequality (|30ll is that, just as 
in thermoelasticity, the static heat flux vanishes: 

Theorem 3.2 The heat flux q vanishes for all termal equilibrium states 
(F, E*^ 0) G V, that is, 

q(F, ^,E^, 0) = 0. (37) 



3.4 Use of Invariant response functions 

In order to satisfy the principle of material objectivity [1] , [5], the functions 
e and ijj, must be a scalar invariant under rigid rotations of the deformed and 
polarized body. The invariance of ^/^ in a rigid rotation is assured when is 
an arbitrary function of the referential quantities E^m, Wl, Gl, where 

Elm = -{Clm — Slm), Clm = Xk,LXk,M , (38) 
Wl = = W = F^E*^ (39) 



OXl dxP OXl 
dG dG dxP 



G = F^g . (40) 



OXl dxP OXl 
Hence we assume that 



^ = V^(E, ^,W, G). (41) 

Next we calculate the time derivatives in equation ( |3T1) by using ip in 
place of ip- We find 

dip p [ OErs dip OWr dip dGn i d dx^ 

W ^ ^OErs d{dx^ /OXk)^ dWRd{dx^ /OXk)^ dGnd^dx^ /dXKyJtdXK 

(42) 
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Now, by dMD-dinD, 



dip 



dE 



RS 



d dx^ 



dijj 1 / dx^ \ dx'^ 



dEns didxyOXK) dt OXk OErs 2 V dXs dX^ ^ OXk 



1 ^ dip dx^ ^ dip dx^ \ dx"^ 



dip dx^ dx"^ 



f— ] 

2ydEKsdXs dEuKdXn^dXK dEuKdX^dXK 



F,(43) 



dip dWR d 9x* 



dip 



dWR d{dxydXK) dt dXx dWn 
and, similarly with the latter, 

dip dGji d dx^ dip 



SkrE'' 



d dx^ 



dip 



dGn d{dx^/dXK) dt dXK dGn 
hence 



' dtdXK ^dW 



d dx^ / dip 

^""""^'jtdx^ = y-dG 



)E 



M 



■F, (44) 



(45) 



moreover, 



dip 
dF 

dip 
de: 



dip J. dip 
dE dW 



.M , dip 



dG 



dip dOn dip . . 

9 = -^;7^^r9 



dGn dg 

dip 



dG 



R 



dE^^ 



■E 



M 



y^dw)-^ ■ 



(46) 



(47) 



(45 



By recollecting the equalities above we can rewrite the dissipation in- 
equality (fT7|) as 



P 



r/^^pT dip 
dE dW 



pAf , dip 



g -F 



-T • Vv 



dip ■ dip \ • ,j / dip 



de 
1 



dG 



q • g + pTT • E*'^ < . 



Now we apply the method of Colemann-NoU in Remark 13. 11 By the 
arbitrariness of g we have dip/dG = O. Remembering ( l33l) . ( IMI) Eq. ( H9l) 
reduces to 



P 



fdip T, dip 



dW 



E'' - p-^F-'- 



+ F 



- P 

dip 
dW 



77 -E 



AI 



(50) 



By the arbitrariness of the time derivatives F, 6', E*^ and by substituting 
the constitutive relation (1^ we find 



iq-g<0. (54) 

We have proved the version of Theorem 13.11 that employes the objective 
energy response function ip. 

Theorem 3.3 The Dissipation Principle is satisfied if and only if the fol- 
lowing conditions hold: 

(i) the objective free energy response function ^, W, G) is indepen- 
dent of the temperature gradient G and determines the entropy, the Cauchy 
stress tensor, and the polarization vector per unit mass through the relations 
(EIP-(E3^; 

(a) the reduced dissipation inequality [54\ ) is satisfied. 

We point out that Equalities fl^ . ffTUj) and fl3^ yield the following ex- 
pression for the Cauchy stress: 

T = pf||f^ - P ® E^^ ; (55) 
hence for the antisymmetric portion of r we obtain the expression 

= ^ (e*^ ® P - P ® E^'^) , (56) 
that coincides with (3.24) of [2]. 



3.5 Internal Dissipation and Entropy Equality 

The local internal dissipation 6o in a thermoelastic body is defined by [Sj 
(p.ll2) 

5o = Of] — ir divc\) (57) 

P 

Then it is proved that (5o = along every local thermoelastic process. Here, 
with regard to thermo-electroelasticity, we define the internal dissipation inst 
by (!57j) and hence we extend the above theorem by the following 
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Theorem 3.4 Along any local process of B we have 

5o = 0. (58) 

Proof . 

By inserting the energy equation (fT^ in the definition (15 7p we obtain the 
equality 

S^ = 9fj--(pe-T-VY-E^^ ■ pir) (59) 
P 



Now by ( 1161) we find 

er] = -ip + e-eT]-E^^ -TT- E^^ ■ n (60) 



and by replacing the latter into (|59|1 we have 

(5o = ( - + £ - - E*^ ■ TT - E*^ ■ tt) - - T • Vv - E*^^ ■ P7t) , (61) 
thus 

6o = ( -07] -E^^ -tt) + -T . (62) 
V ^ p 

Now by (I3l])i 

5„= (^-d^'^.p-de^-e-dEM^-E^^ -d^i^-g-9r]-E^^ -TT^+^T-Vv (63) 

and the constitutive restrictions fl?r|) - fl2^ together with flH^ yield fl3Hl) . ^ 

In a thermoelastic body any thermoelastic process is locally reversible, in 
the sense that the following entropy equality holds 

r divci 

PV = P-Q Q- ■ (64) 

The last theorem and (l57j) yield the same result in thermo-electroelasticity 
too. 

Theorem 3.5 Along any local process of B the entropy equality holds. 
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4 Referential description 



We can rewrite the constitutive relations (|T8l) -(l22l) in material form by using 
the first Piola-Kirchhoff stress tensor S(X, t), that is related with the Cauchy 
stress by 

S = JF-V (65) 
and by using the well known equalities ([5]), (E])-® and 

DivS = JdivT , (66) 

Now the process class IP{B) of B of Section [21 contatining the processes 
(lH), must be substituted with 1Pr{B), that is the set of ordered 10— tuples 
of functions on B x iR 

= (x(.), e{.), ^(.), e(.), ^(0, S(.), IP(.), Q(-), b(.), r(.)) G FniB) 

(67) 

defined with respect to B , satisfying the material versions of the balance 
laws of linear momentum, moment of momentum, energy, the entropy in- 
equahty,and the field equations of electrostatics. 

4.1 Local Balance Laws in Material Form 

Under suitable assumptions of regularity, and using (jl]), the usual integral 
forms of the balance laws of linear momentum, moment of momentum, en- 
ergy, the field equations of electrostatics, and the entropy inequality are 
equivalent to the referential field equations 

p^v = DivS + p/jb , (68) 

Pre = S-F- DivQ + W ■ IP + p^r , (69) 

W = -V^^ (= -F^V^(p) , DivA = , (70) 

PbV > pn{r/e) - Dzv{Q/9) . (71) 
Let ip = ipO be the specific free energy per unit mass defined by 

^ = e-e7]~W -U. (72) 
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Then f lU^ and flTTl) yield the dissipation inequality 

Pi?(V' + ^^)-S-F + ^Q-G + IP-W < 0, (73) 

where G = Grad6(K, t) is the referential temperature gradient. 

Remark 4.1 Note that the free- energy function ip defined here by f/^ ) co- 
incides with the analogous function x defined in (4-2) of on page 596. In 
fact, 

W • n = W • W/pn = {FV) ■ {JF''P)/pn = E^^ • P{J/pr) 
Hence, by pr = J p, we have 

W-n = E^^-P/p = E^'^-7r, 

where tt is the spatial polarization vector. 

4.2 Referential Constitutive Assumptions 

Let Vr be the open, simply connected domain consisting of 4— tuples (F, 6, W, G) 
such that (F, 9, W, G) G V- hence if (F, 9, W, G) G Vr, then (F, 9, W, 0) G 
Vr. 

Next we use a free energy function of the form 

= ^(F, 9, W, G) . (74) 

Assumption 4.1 For every p G 1Pr{B) the specific free energy ip^X., t), 
the specific entropy ?7(X, t), the first Piola-Kirchhoff stress tensor S(X, t), 
the specific polarization vector IP(X, t), and the heat flux Q(X, t) are given 
by continuously differentiable functions on T>r such that 





= ^(F, ^,W, G), 


(75) 


V 


= r)(F, ^,W, G), 


(76) 


s 


= S(F, ^,W, G), 


(77) 


IP 


= IP(F, 9,W, G), 


(78) 


Q 


= H(F, ^,W, G). 


(79) 


Of course, once Pr{.), IP(. 


,), ^/;(.) and i7(.) are 


known, then equality fl72|) 



gives the continuously differentiable function e{.) determining £:(X, t) such 
that 

e = e{F,9,W,G). (80) 

Also note that the dependence upon X is not written for convenience, but 
it is implicit and understood when the body is not materially homogeneous. 
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4.3 Coleman-Noll Method and Thermodynamic Re- 
strictions 



Given any motion x(X, t), temperature field ^(X, t) and electric poten- 
tial field v^(X, t), the constitutive equations (!73I) - (17H1) determine e(X, t), 
r/(X, t), S(X, t), IP(X, t), and the local laws (EHD and (El) determine 
b(X, t) and r(X, t). Hence for any given motion, temperature field and 
electric potential fielf, a unique process p is constructed. 

The method of Coleman-Noll [3] is based on the postulate that every 
process p so constructed belongs to the process class IP{B) of B, that is, on 
the assumption that the constitutive assumptions ( I75l) -( !79|) are compatible 
with thermodynamics, in the sense of the following 

Dissipation Principle For any given motion, temperature field and 
electric potential field, the process p constructed from the constitutive equa- 
tions ([y3y-([yP|j belongs to the process class IP{B) of B. Therefore the con- 
stitutive functions [75\)-l^{7m are compatible with the second law of thermody- 
namics in the sense that they satisfy the dissipation inequality ( [77p . 

Theorem 4.1 The Dissipation Principle is satisfied if and only if the fol- 
lowing conditions hold: 

(i) the free energy response function ip{F, ^, W, G) is independent of the 
temperature gradient G and determines the entropy, the first Piola-Kirchhoff 
stress, and the polarization vector through the relations 

r)(F, e,W) = -9,^(F, ^,W), (81) 

S(F, e, W, Q) = Pk9f7A(F, W, Q) (82) 

fl(F, e, W, Q) = -5wV^(F, e, W, Q) (83) 

(ii) the reduced dissipation inequality 

Q • G < (84) 

is satisfied. 

Proof . By he chain rule we have 

^ = d^^ ■¥ + deiJ - e + dw-^ ■ W + Og^ ■ G . (85) 

Thus by substituting this equation together with the constitutive equa- 
tions (!75l) -( |79l) into the dissipation inequality (|73ll gives 
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{pnd^i; - S) ■ F + {pRdei^ + fi)e + {pRchfi^ + IP) ■ W 



(86) 
(87) 



Now we use Coleman-Mizel method j3]: by Remark l3. 11 translated in ref- 
erential form, we have that F, ^, W and G can be assigned arbitrary values 
independently from the other variables and the theorem is easily proved just 
as in the proof of Theorem 13.11 <^ 

A consequence of the reduced dissipation inequality is that, just as 
in thermoelasticity, the static heat flux vanishes: 

Theorem 4.2 The heat flux Q vanishes for all termal equilibrium states 
(F, e, E*^ 0) e V, that is, 

Q(F, ^,W, 0) = 0. (88) 
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